The Implications of Discontinuities for Testing Theories of Turbulence in the Solar 

Wind 
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In-situ observations of magnetic field fluctuations in the solar wind show a broad continuum in the 
power spectral density (PSD) with a power-law range of scaling often identified as an inertial range 
of magnetohydrodynamic turbulence. Ifowever, both turbulence and discontinuities are present in 
the solar wind on these inertial range of scales. We identify and remove these discontinuities using 
a method which for the first time does not impose a characteristic scale on the resultant time-series. 
The PSD of vector field fluctuations obtained from at-a point observations is a tensor that can in 
principle be anisotropic with scaling exponents that depend on background field and flow direction. 
This provides a key test of theories of turbulence. We find that the removal of discontinuities from the 
observed time-series can significantly alter the PSD trace anisotropy. It becomes quasi-isotropic, in 
that the observed exponent does not vary with the background field angle once the discontinuities are 
removed. This is consistent with the predictions of the Iroshnikov-Kraichnan model of turbulence. 
As a consistency check we construct a surrogate time-series from the observations that is composed 
solely of discontinuities. The surrogate provides an estimate of the PSD due solely to discontinuities 
and this provides the effective noise-floor produced by discontinuities for all scales greater than a 
few ion-cyclotron scales. 
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Turbulence is ubiquitous in astrophysics and may 
be important in cosmic ray scattering [ll-Q. In-situ 
satellites provide at-a-point observations of heliospheric 
plasma parameters over a broad range of temporal scales. 
Solar wind observations display fluctuations from plasma 
kinetic to solar rotation scales Q and exhibit ranges of 
power-law scaling. Evidence for Magnetohydrodynamic 
(MHD) turbulence includes: observed piece-wise power- 
law scaling in both the components and the trace of 
the magnetic field power spectral density (PSD) [H-Q, 
non-Gaussian probability distribution functions (PDF) 
of the fluctuations that display multi-fractal character- 
istics [H, H, and a large estimated effective magnetic 
Reynolds number l^, 11|. There has as a consequence 
been considerable interest in using the solar wind to test 
the predictions of MHD theories of turbulence. 

In studies of solar wind turbulence, the models of 
Iroshnikov-Kraichnan (hereafter IK) [l^, and Goldre- 
ich and Sridhar (hereafter GS) [13] have merited partic- 
ular attention (e.g. (gI [l5l - [25j ). The isotropic IK model 
of MHD turbulence predicts E{k) ^ k^^/"^ , whilst the 
anisotropic model of GS predicts E{kj_) ~ k^ ' and 
E{k\\) fcjj"^, where E{k) is the energy spectral density 
as a function of Fourier wavenumber, k. The compo- 
nents k± and fcy refer to the orientation of the wavevec- 
tor w.r.t to a background magnetic field direction, which 
is often interpreted to be local and scale dependent (e.g. 



curately. Distinguishing between the predictions of these 
two theories from data is highly topical (e.g. [1,@,[22, 26- 

m). 

In the solar wind there are discontinuities in addition 
to the turbulence. Quantifying the effect of discontinu- 
ities on statistical measures used in the framework of 

ii-IHi. 
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[d, Ullii, [iai). Thus, in principle in-situ observations 

of the solar wind offer a " laboratory" to distinguish be- 
tween the distinct predictions of IK and GS providing it 
is possible to identify the background field direction ac- 



turbulence is a challenge to data analysis 
portantly, discontinuities in a time-series can result in 
power-law scaling with exponents in the range predicted 
by turbulence theories [sj] . Accurate identification of the 
background magnetic field direction is critical to testing 
theories of anisotropic MHD turbulence. However, this 
can be compromised by the presence of discontinuities. 
In this Letter we employ a new method to identify and 
remove fluctuations affected by discontinuities. Disconti- 
nuities are usually identified by a change in the magnetic 
field magnitude and/or direction between two observa- 
tions that exceeds a given threshold 31, Hi- Uniquely, 
our method does not assume a single characteristic phys- 
ical scale for this threshold and hence for discontinuities 
in the solar wind, and so, does not impose a single scale 
upon the observed time-series. This is particularly im- 
portant when investigating turbulence, which by defini- 
tion has no characteristic scale. We find that, once these 
discontinuities are removed, the remaining signal is con- 
sistent with IK but inconsistent with GS, specifically, the 
trace PSD exponent is isotropic with an exponent close 
to -3/2 for both k±_ and fc|| sampling directions. 

We focus on observations of fast, quiet solar wind us- 
ing the Ulysses satellite [day 91-146, 1995] whilst at a 
heliographic latitude of 21° — 58° and a radial distance 
1.36 — 1.58 AU. The interval is recorded at a cadence of 
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Figure 1: P^f^j^ (red), Pf^p (blue) and Pib\,r (black). C, R, F 
denote component, raw and discontinuity filtered respectively 
(See text for details). The spacecraft frequency is on the bot- 
tom X-axis with the spacecraft time-scale along the top x-axis 
with blue and red type for minutes and seconds respectively. 



1 s in an average solar wind flow of 756 km/s. To cor- 
rect for data gaps and inhoniogeneous sampling the data 
is linearly interpolated to a cadence of 2 s. For periods 
of missing data longer than 10 s the resulting region of 
interpolated data is excluded from subsequent analysis, 
thus limiting the effect of data gaps and interpolation. 

Any study of anisotropic plasma turbulence requires 
a definition of the background magnetic field direction. 
The background field, Bq, is defined by an average over 
some time-scale r. Previous studies have considered a 
single time-scale, for example, that of the outer-scale of 
the turbulence (e.g. [s^l)- Alternatively, a local back- 
ground field is defined as a scale dependent estimate over 
a Twhich now varies, where r relates the scale of aver- 
aging to the scale at which fluctuations are calculated, 
either using GSF or wavelets (e.g. d, H^). 

The method we will use is based on an assumption 
that the magnitude of every vector in the time-series 
may be written as = (Bo + + ^l: where Bq is 
the background field and || are some small addition 
to the background field. Taylor expansion of this ex- 
pression to lowest order leads to (5|B| = |Bi| — IB2I ^ 
^0,1 + ~ ^0,2 — e||,2- Thus, if the background field is 
unchanged between two observations Bq i = i3o,2, then. 
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SB, 



where SB^^ — 
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Therefore, in the 



absence of discontinuities in the observed time-series the 
fluctuation power of the modulus of the magnetic field, 
P\B\j and the power in the parallel component of the 
fluctuations w.r.t the background field, , will be ap- 
proximately equal. This suggests an important criterion 
for determining whether discontinuities have a significant 
impact on the statistical properties inferred from the ob- 



servations. Note, unlike , P|b[ is independent of any 
averaging scheme used to estimate the background field. 
Therefore, the difference between P^ and P|b| may be 
regarded as a measure of the quality of the background 
field estimate. 

Here, we use a local background magnetic field esti- 
mate, obtained for each scale r by averaging the observed 
time-series with a Gaussian window of width 4r, where 
the standard deviation of the Gaussian window deter- 
mines the spectral characteristics in Fourier space and 
we may define / = l/r [iij. We define the method for 
detecting discontinuities as follows: a discontinuity is ob- 
served at time, t, if: 



|<5|B(t,/)|-<5i?|| {tJ)\>Tor 



(1) 



where Tpand /3 are values to be found from the data that 
determine a scale dependent threshold for the discontinu- 
ities. We have explored a broad range of Tg and /3. Here, 
we show results using Tq — 0.022 and (3 = 0.61, which 
have been selected due to best agreement between P^'p 
and P\B\,Rj where we refer to the original observations 
as "raw"(R) and observations filtered by equation ^ as 
"discontinuity-filtered" (F). 

We first show that our procedure removes discontinu- 
ities in a consistent fashion. In Fig. [T]we plot the PSD 
of P^'j^ (red), P^p (blue) and P\b\.r (black). We can see 
that, > P\B\,R (i-e- > S\B\) and the spectral 
indices are distinct for time-scales t < 30 min. However, 
when we consider the discontinuity-filtered case we can 
P\B\.R within experimental uncertainty. 



see that P^p 



Agreement between P^p and P\b\,r justifies our choice 
of To and (3. Thus, we have removed the effect of discon- 
tinuities from the observed time-series of the components 
of the magnetic field in the solar wind according to our 
definition ([T]). 

When observations are made at a single point in 
the solar wind it is only possible to infer the re- 
duced one-dimensional energy spectra. Thus, the ob- 
served spectral tensor may be expressed as Pij (/, 9) ^ 
J d3kPy (k)(5(27r/ - k • V,„), where P^jif,e) is the ob- 
served spectral tensor, Pij(k) is the three-dimensional 
spectral tensor, Ysw is the solar wind velocity, k is the 
wave vector, / is the frequency in the spacecraft frame 
and 9 is the angle between the sampling direction and the 
local background magnetic field 35| . The trace power is 
then Pj{f) = l^iPii{f,9). In the case where 9 is subdi- 
vided such that 9 = 80° -90° ~ 9± and 6* = 0° - 10° ~ 6i|| 
the corresponding trace power may be regarded as an es- 
timate for the spectra of k± and fc|| respectively. For 
IK it is expected that n (/) - /-3/2 and for GS it 
is expected that Pj(/) - /-^/^ ^^^^ pT(^j^ _ j-2 

the solar wind, P^ has been found to be anisotropic, 
in that the power and spectral index vary with 9 (e.g. 
[23, 26 1). The scaling is not always consistent with 
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—5/3 and —2, but generally there is agreement that the 
spectral indices of — '^1/3 [26[. These observa- 
tions are interpreted to support the GS model (See e.g. 
Refs dm III). 

We now analyse the raw time-series alongside the 
discontinuity-filtered time-series to determine how dis- 
continuities affect the measured trace power of the so- 
lar wind when sampling parallel, , and perpendicular, 
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j_ , to the local magnetic field. We plot the PSD of 
the trace power in Fig. [2] In this figure the squares 
and circles represent the case of Pj^ and P^ . The left 
and right panels show the raw and discontinuity-filtered 
cases respectively. The solid black lines show fits with 
PSD exponent —1.9 and —1.5 where appropriate. 

Both P|j^j^ p (RH plot) are at lower power than their 
raw counterparts (LH plot), especially at large scales. 
Remarkably, the effect of removing discontinuities is to 
transform the exponents of the PSD from anisotropic (LH 
plot) to isotropic (RH plot). Once discontinuities arc re- 
moved, the spectral exponents are approximately equal 
for both P^ and Pf , with a value of -I.5±0.05. This re- 
sult is dramatically different to that suggested by the raw 
observations that show that the spectral exponent of P'^ 
is greater than Pj, with a difference of ~ 1/3 [26|. Also, 
in this interval there is no data above t ^ I hr, imply- 
ing an upper limit to the time between discontinuities, 
irrespective of sampling direction w.r.t the background 
magnetic field. 

Our results show that the presence of discontinuities 
can significantly affect the anisotropy seen in both the 
power and scaling exponents of the PSD. Thus, any quan- 
titative comparison between observations and theories of 
turbulence must address the issue of discontinuities. Fur- 
thermore, using our procedure we obtain a result that 
supports IK scaling, similar to Ref. [25]. Intriguingly, we 
also find good agreement with a recently propose phe- 
nomenological theory of incompressible MHD turbulence 
28l 36 1 . However, it is not possible to completely confirm 
this new phenomenology from this study as the parallel 
and perpendicular dissipation scale must be measured. 
This is not possible with the available cadence of the in- 
struments used in this study. 

Finally, we will demonstrate that a simple surrogate 
time-series generated entirely from discontinuities can ac- 
count for some, but not all, of the power in the in-situ 
observations of fluctuations in the solar wind. We will 
construct a surrogate dataset from the Ulysses fast solar 
wind observations. 

The surrogate data is a Continuous- Time Marko- 
vian jump process, where the angular state space of 
the observations is preserved in the surrogate time se- 
ries. The observations, in RTN coordinates ( Bt = 
[i?fl(i), i?T(i), 5Af(i)]), are projected onto spherical co- 
ordinates (Bi = [|P|f , 6*4, 0t]), we may then define the 
joint probability, Pe,^, between angles 9 and (j). We will 
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Figure 2: Trace power Pj^ and P"^ shown by squares and 
circles. Left panel: Raw case. Right panel: Discontinuity- 
filtered case (See text for details). Examples of power index 
—3/2 and —1.9 are shown where appropriate with solid black 
lines. The spacecraft rest frame frequency is on the bottom 
X-axis with corresponding time-scale along the top x-axis with 
blue and red type for minutes and seconds respectively. 



construct our surrogate time-series to have the same Pg ^ 
as the original observations. 

The waiting time between jumps, Tw, is an i.i.d ran- 
dom variable. The waiting time distribution of the 
solar wind has been measured by numerous authors 
II H S m, 113, [si]. Reported distributions are ex- 
power-law and log-normal. A close approx- 



ponential, 

imation to all of these probability distributions is the 
Weibull (Stretched-exponential) distribution [s^, such 
that Twlt\\,b) = !^{!-f-^eyi]y{-{{f). Thus, we use 
the Weibull distribution as our waiting time distribution. 
The values for our distribution are chosen by producing 
the best fit between the PSD of the observations and the 
surrogate in the time scales, r ^ 1 hour to r ~ 7 days. 
This is achieved by A = 75 seconds and b = 0.4. With 
these chosen parameters the mean waiting time is ap- 
proximately 7 min. As our surrogate time series is a jump 
process the value of our surrogate time series at a specific 
time, Xt, remains constant until a jump occurs. When 
a jump occurs a random set from the observations are 
selected, Bj, where j is a randomly selected index with 
uniform probability for all data points, which effectively 
randomises the phase of the discontinuity fluctuations, 
whilst preserving the observational Pe,<^. Our surrogate 
time series may be written as 



lMo,0(x;),0(a:;)], 



iit = 

if < > and t - s < Tw 
iit>0 andt- s^Tw 
(2) 

where 1/A is the jump rate, t is the time of the surro- 
gate data, s is the time of the latest jump in the time 
series, 9 (x^) and (j> (a;^) are subsets from the same time 
index of the observations. Mq is the constant magni- 
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Figure 3: Comparison between the observations (red) and 
surrogate data of discontinuities (blue) (See text for details). 
For clarity only the R component is shown. Spacecraft fre- 
quency is along the bottom, whilst spacecraft time-scale is 
along the top in minutes (blue) and second (red). 



tude of the vector, chosen such that there is agreement 
between the PSD of the observations and the surrogate 
date in the range t ~ 1 ~ 7 days when the surrogate has 
been projected back to RTN coordinates (these scales are 
not shown here to improve clarity of the inertial range 
scales). The surrogate data is constructed in spherical 
coordinates, but for the purposes of comparison is plot- 
ted in RTN coordinates. 

We plot the PSD of the surrogate time-series (blue) 
alongside the observations (red) in Fig. [3l We can see 
that the observations are closely reproduced by the sur- 
rogate at scales r > 5 minutes. However, where there is 
no agreement within uncertainty there must also be an- 
other process with different autocorrelation properties in 
the observations distinct from that of the time-series of 
discontinuities, e.g. turbulence. This can be seen in Fig. 
|31 Importantly, we do not claim that only discontinuities 
are present on time scales r > 5 minutes, only that dis- 
continuities are increasingly the dominant contribution 
in the observations on these larger scales. 

We have presented evidence for more than one pro- 
cess in the solar wind inertial range. Discontinuities 
show increasing PSD dominance as the scales increase 
but can not account for all the smaller scale fluctuations 
of the inertial range. Thus, there must be at least two 
physical processes present in the inertial range, i.e. dis- 
continuities and turbulence. The challenge is to sepa- 
rate the PSD contribution from turbulence and discon- 
tinuities. We achieve a separation by using a thresh- 
old that does not introduce a characteristic scale to the 
time-series under investigation. We measure the differ- 



ence between two estimates of the parallel component 
((5J5|| and S\B\) to exclude data where the local field is 
poorly defined/estimated. We show that our method re- 
moves discontinuities, as P^p — P\b\,r- The trace power 
is shown to to scale in a manner that is independent 
of Obv, this is consistent with the IK and more recent 
quasi-isotropic cascade phenomenology 28, 3^. Our re- 
sult does not support the GS model in the solar wind, but 
suggests previous agreement may be due to a poorly de- 
fined local field due to discontinuities in the time-series. 
We purpose a new model for discontinuities in the solar 
wind that demonstrates the contribution to all scales in 
the coordinate system of the observations. We find that 
discontinuities must be considered when using PSD ex- 
ponents to test theories of turbulence in the solar wind. 

The authors acknowledge the Ulysses instrument 
teams for providing magnetometer data. This work was 
supported by the UK STFC. 
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